Low friction and rotational dynamics of crystalline flakes in solid 
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Abstract. - Solids at incommensurate contact display low- friction, 'superlubric', sliding. For 
graphene flakes on a graphite surface, superlubric sliding is only temporary due to rotation of the 
flakes from incommensurate to commensurate contact with the substrate. We examine this rota- 
tional channel of friction in a prototype geometry of meso- and macroscopic solid lubrication. By 
molecular dynamics simulations and theoretical arguments we find that two surfaces lubricated by 
mobile, rotating graphene flakes exhibit stable superlubric sliding as for ideally incommensurate 
contacts also when they are covered by randomly oriented pinned graphene patches. For com- 
mensurate surfaces, we find a low friction state at low temperature where incommensurate states 
are not destroyed by thermal fluctuations. 



Superlubricity between incommensurate surfaces pro- 
vides a desired low-friction state essential for the function 

' of small-scale machines and function of solid lubricants. 
Vanishing static friction has been first predicted by Aubry 
and Peyrard [1,2 for infinite lattices on an incommensu- 

. rate periodic potential. Later, Shinjo and Hirano [3] pre- 
dicted that for infinite incommensurate contacts also the 

. kinetic friction would vanish and called this effect superlu- 
bricity. Extremely low friction has been observed for small 
contacts in AFM experiments by Dienwiebel et al. [4] for 
the sliding of small graphite flakes on a single surface of 
graphite. However, the superlubric sliding is experimen- 
tally found to be only temporary, going over to stick-slip 
behavior after several scans along the surface [5]. Direct 
molecular dynamics simulations [5] as well as a theoret- 
ical analysis of the nonlinear dynamics [6] demonstrated 
that this increase of friction is due to the rotation of the 
flakes to a commensurate contact with the single surface. 
The incommensurate orientations of the flakes are stable 
but not robust against thermal fluctuations particularly 
if the center of mass of a flake travels on top of the sub- 
strate atoms [6]. The rotation of nano-scale crystals such 
as graphene flakes is the simplest nontrivial issue in lubri- 
cation with lamellar solids such as graphite and M0S2. 




Fig. 1: (Color online) The system describes two infinite 
plates with graphene flakes embedded in between. The plate 
surfaces are covered by single-domain or randomly oriented 
multi-domain graphite layers, as shown for the bottom 
substrate. The top plate is pulled by a spring, the tension 
in which gives a measure of friction. The pulling angle 
relative orientation of lattices on the two plates f3, and flake 
orientation 6 are indicated. 



In this Letter, we propose a model system that can be 
used to investigate solid lubrication at both nano and 
macro scales, and which exhibits stable super- low fric- 
tion. The system, shown in Fig. [TJ consists of two infi- 
nite plates covered either by single-domain or randomly 
oriented multi-domain (graphite) layers with (graphene) 
flakes embedded between them that act as a solid lubri- 
cant. As such, the model contains all key elements of 
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systems lubricated by lamellar solids, namely a double in- 
terface, multiple sliding objects, and rotational dynamics, 
which are missing in previous models of friction of a single 
flake on a single surface [5j[6] and in the models of lubri- 
cation with atoms and molecules [7H9]. Conversely, we 
simplify the interaction compared to more precise quan- 
tum chemistry calculations such as those of Ref. [10 j. The 
flakes could be inserted intentionally or result from wear 
of lamellar solids, as was found in [TTHT3] . We show 
that macroscopic low-friction sliding, expected for single- 
domain incommensurate contacts, also occurs for disor- 
dered plates covered with randomly oriented domains once 
lubricating flakes are introduced between them. In fact, 
although the lubricant flakes readily get into commensu- 
rate contact with either of the plates, stick-slip behavior 
is prevented because the flakes cannot be commensurate 
with both plates at the same time, contrary to adsorbed 
flexible chain molecules ('third bodies') [7-9 . Besides this, 
by examining the stability of the superlubric sliding, we 
find a non trivial mechanism yielding low friction also for 
commensurate surfaces at low temperature. 

Although our findings are based on a simplified model, 
considering superlubricity mediated by the flakes gives a 
perspective to exploit a genuinely atomistic mechanism at 
larger length scales. Multi-domain graphene layers can 
be already obtained on macroscopic scales [14 and sev- 
eral other synthesis routes are being currently pursued, 
making this approach suitable for realistic applications. 
Disordered multi-domain surface structures can also be 
formed spontaneously due to capillary forces which glue 
randomly oriented flakes to the sliding surfaces, resulting 
in a random carpet of flakes rather than a single infinite 
domain. This effect prevents rotation of flakes and en- 
sures persistent incommensurability. Our model and the 
results have important implications for understanding the 
macroscopic properties of graphite and other solid lamel- 
lar lubricants, which are undoubtedly the most common 
solid lubricants [15,16 . Transmission electron microscopy 
(TEM) observations on MoS 2 and WS 2 [EJUS] have in- 
dicated that in macroscopic sliding contacts of lamellar 
solids rotated flakes are created. 

We consider n rigid hexagonal graphene flakes formed 
by N atoms each, which are free to move and rotate in 
the plane between the two plates. The typical coherence 
length of graphene, which is about 1/im, is much larger 
than the size of a flake, and therefore elastic deformations 
can be neglected [T71IT9] . The coordinates of a flake i 
are defined by the center-of-mass position and the ori- 
entation 4>i with respect to the x-axis. In the numerical 
simulations, we have taken the total number of embedded 
atoms nN = 3456, which corresponds to sufficiently large 
number of flakes. Further increase of their number does 
not affect essentially the friction response. 

In order to mimic the commonly used experimental fric- 
tional setup, we consider a model where one of the surfaces 
is stationary, while the other is pulled by a spring coupled 
to a support moving at constant velocity V. We consider 



friction between atomically smooth surfaces, which is the 
typical configuration in surface force apparatus (SFA) ex- 
periments where the scale of the roughness is larger than 
the size of the flakes so that it does not influence the trans- 
lation and rotation motion of the flakes, and the mecha- 
nism described in this Letter contributes significantly to 
the friction. 

Each flake atom is subjected to a potential due to the 
interaction with each of the two surfaces that consist either 
of single domains or of patches with random orientation. 
The potential function has hexagonal symmetry 

rr I \ Inx' 2th/ 4th/ \ 

(J (x, y) = —Uq I 2 cos — — cos— h cos — - — 1 , (1) 

\ Ai A2 A 2 J 

where (V, y') = (x cos /3 + y sin /3, —x sin f3 + y cos 0) is the 
position of a flake atom on a graphite domain with orien- 
tation f3. The parameters Ai = 0.246 nm, A2 = \/3Ai, and 
Uq = 2.26 meV produce the structure of graphite with a 
corrugation amplitude of 10.17 meV. 

The dynamics of the system is modelled using Langevin 
equations which include random forces rji(t) and fji(i) and 
random torques Q(t) and (i(t) : to describe the effects of 
thermal fluctuations. Tildes are used to indicate quanti- 
ties related to the moving plate. The equation of motion 
for the position R of the moving plate is coupled to those 
of the flakes as 

n 

MR = ^ [-Fflake(r* - R, 0<) + jNm(ii - R) 
1=1 

n 

-«(R-Vt) (2) 

i=l 

NmYi = FflakeOi, (j>i) + FflakeOi - R, 0») 

- jNmvi - jNm(ri — R) + rji(t) + fji(t) , (3) 

I<j>i = Tflakefo, fa) + T flake (r i - R, (t>i) 

+ + , (4) 

where M is the mass of the driven plate and m is the mass 
of a carbon atom, / is the moment of inertia of a flake, 
Fflake( r , <j>) [Ffl a ke( r 7 <j>)\ is the force acting on a flake at 
position r relative to the plate and orientation <p due to 
interaction of all atoms with the surface of the stationary 
(moving) plate, and T fla ke(r, <j>) [T flake (r, 0)] is the corre- 
sponding torque. Similar calculations can be done also for 
multilayer flakes. What is essential is the structure of the 
layers which are in contact with the substrates. 

The lateral force measured in friction force experiments 
F = k;(R — Vt) is the force acting on the spring when the 
support is moving at constant velocity V and its time aver- 
age is equal to the friction force Ffriction = | (F) | . Molecular 
dynamics simulations of the system have been performed 
using a fourth-order Runge-Kutta integration algorithm. 
The precision of the algorithm has been tested by energy- 
conservation at T = without driving. We point out 
that if energy conservation in this case is not obeyed suf- 
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Fig. 2: Lateral force, normalized to Fk = ftAi, as a function of 
the support displacement with N = 24. The dashed line is the 
lateral force for two commensurate surfaces each with a single 
infinite domain. The dotted line is the force for two plates 
with infinite domains at a relative incommensurate orientation 
of P — 30° . The solid line is for two surfaces covered by random 
domains with a domain size of 2.46 nm. For two commensu- 
rate surfaces, the friction shows distinctive stick-slip behavior, 
whereas for the surfaces with domains and for incommensurate 
ideal surfaces the forces do not display stick-slip behavior and 
are of similar order of magnitude. 



ficiently, the stability of the incommensurate orbits is en- 
hanced, artificially reducing the friction. With the fourth- 
order Runge-Kutta algorithm energy-conservation is en- 
sured on the time scale of our simulations with a time 
step of 0.02 ps. 

The moving plate of mass M = nNm, with rn = 
1.99 x 10 -26 kg the mass of a carbon atom, is pulled by a 
spring with spring constant k = uNhcq = 4.15 nN/nm cou- 
pled to a support moving at constant velocity V = 12 m/s 
at an angle of 6 = 70° with respect to the x axis. The 
domains are 10Ai in diameter and the borders between 
domains are smoothed out by a sinusoidal function. We 
have checked that a reduction in support velocity by up 
to three orders of magnitude lower does not qualitatively 
affect the results. We have also verified that the pulling di- 
rection and domain size affect the results only marginally. 

In Fig. [2j we show the instantaneous friction forces cal- 
culated for the system where the flakes are embedded be- 
tween fully commensurate plates (dashed line), incommen- 
surate plates (dotted line) and plates with disordered do- 
mains (full line). One can see that only for fully commen- 
surate plates there is stick-slip behavior resulting in high 
friction, while for incommensurate and disordered plates, 
the friction forces are comparable and drastically lower. 
This is due to the fact that the flakes cannot lock into 
commensurate contacts with both plates, not only when 
the plates are incommensurate, but also in the presence 
of patches. Thus the lubrication properties of flakes differ 
significantly from those of short flexible chains of particles 
studied in Refs. [7H9], where it is shown that, even for in- 
commensurate plates, short chains can lock to both plates 
and increase the friction. This feature, together with the 
fact that wear in lamellar systems leads to creation of lu- 
bricant flakes may explain why layered materials p~5|fT6] 
are excellent lubricants. 
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Fig. 3: Orientations fa of 24- atom flakes (left) and their dis- 
tribution functions (right) at T = 77 K, for two commensurate 
surfaces [(a) and (b)], two incommensurate surfaces at f3 = 30° 
[(c) and (d)], and for two multi-domain surfaces [(e) and (f)]. 
Note in (a), (b), (c), and (d) that the flakes have orientations 
commensurate to either of the two graphite surfaces. In the 
case of (c) and (d), the friction is low, because the flakes are 
not commensurate to both surfaces at the same time. This 
picture is preserved for multi-domain surfaces. 



Under humid conditions, the multi-domain surface 
structures can be formed spontaneously due to the cap- 
illary forces which fix randomly oriented flakes at the slid- 
ing surfaces, while in vacuum conditions the formation of 
such structures is less likely. In addition, the water vapor 
trapped between the surfaces may saturate the bonds on 
graphite, thus preventing covalent bonding between the 
graphite surfaces. 

These may be the reasons why graphite is such a bad 
lubricant in vacuum, and needs the humidity of air to 
perform well. 

Since the main mechanism of friction is determined by 
the rotation of the flakes from incommensurate (superlu- 
bric) to commensurate (stick-slip) orientation, the ensem- 
ble of angles fa of the flakes provides information on the 
dynamics. We find the surprising result that, in all three 
cases, for most of the time, the flakes are in commensurate 
contact with either of the plates. However, except in the 
case of commensurate plates, the flakes cannot be com- 
mensurate to both plates at the same time. It is precisely 
this frustration that leads to the low value of the fric- 
tion. As shown in Fig. [3j bands at commensurate angles 
are clearly visible for both single-domain commensurate 
[(a) and (b)] and incommensurate [(c) and (d)] contacts, 
whereas the disordered patches [(e) and (f)] ensure that 
all orientations are present. That the latter situation is fa- 
vorable for low friction can be understood from Fig. H^a), 
where we show the dependence of friction on the relative 
orientation of the sliding plate to the stationary one for 
different flake sizes. Only in a narrow range of angles /3 
close to perfect commensurability the friction is high. For 
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Fig. 4: Friction force for various flake sizes as a function of 
the relative angle of the plates at T — 77 K (a), as a function 
of temperature for commensurate plates (b), and as a function 
of temperature for different rotational dynamics of 24-atom 
flakes embedded between commensurate plates (c). In panel 
(a), results for 24-atom flakes between two plates with patches 
are included for comparison and indicated with a solid line at 
^friction / Fk = 0.54. This value is close to the average over all 
angles for single-domain incommensurate plates, Ffriction/^fc = 
0.46 ± 0.05. 



all other angles, the friction is low, except in the case of 
single particle lubricants (N = 1), which cannot create in- 
commensurate contacts [7J[9]. For multi-domain surfaces, 
the relative orientation of the two lattices is random, and 
therefore, the friction will be equal to an average over all 
relative angles, resulting in a large reduction with respect 
to the commensurate contact. 

In Fig. HJb) we show that above a given temperature, 
the friction decreases, as expected. This decrease sets in 
later for larger flakes. Again, patches on both plates lead 
to much lower friction at all temperatures. 

In Fig. ffl^c), we shown that for commensurate plates, 
the temperature dependence of the friction exhibits a 
crossover around T=200 K. Above this temperature, 
the friction drops, following a power-law proportional to 
T _113±0 04 . This crossover can be understood by com- 
parison with the friction calculated by freezing the orien- 
tation of the flakes either in the commensurate orienta- 
tion 4>i = 0, or with a random distribution of angles. At 
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Fig. 5: Role of incommensurate orbits at low temperatures at 
2.6 K and 7 = 0.4/ps. Lateral force (a) and flake orientation 
[(b) and (c)] as a function of support displacement. The distri- 
bution of orientations after long time is shown in (d). Rotations 
from incommensurate to commensurate orientations and vice 
versa are caused by the energy released during a slip. In be- 
tween slips, the incommensurate orientations survive, because 
they are in principle stable and they contribute significantly 
to the total distribution of orientations, and therefore lead to 
lower friction. 



low temperatures, most flakes have <^ = 0, whereas at 
higher temperatures, incommensurate orientations come 
into play and the friction is similar to that of a frozen 
system with flakes with random orientations. 

Lastly, we discuss the unusual reduction of friction 
found for all flake sizes at very low temperatures [Fig. [4(b) 
and Hl[c)] . This effect is due to the existence of stable in- 
commensurate orbits and is shown in more detail in Fig.JU 
These orbits can be present in the initial conditions, but 
they are also activited during slips. For two commensu- 
rate plates a great deal of energy is released into the sys- 
tem during a slip [see Fig.[5ja). This causes some flakes to 
move into the incommensurate orientations [see Fig. [5f b) 
and (c)], which are not robust against thermal fluctua- 
tions. At high temperatures, the incommensurate orbits 
decay rapidly to commensurate ones (<j) — 0, 60, 120°), and 
so their effect on friction is negligible [5 ,6 . At low tem- 
perature, however, the decay rate of the incommensurate 
orbits is so low, that they survive for a long time, signif- 
icantly reducing the friction. The distribution of orienta- 
tions is shown in Fig.[5fd), where one can see peaks at the 
commensurate orientation, as in Fig. [3fb), and at incom- 
mensurate orientations at <j) = ±26, 34, 86, 94, 146, 154°. 
These coincide precisely with the stable incommensurate 
orientations found for N = 24 [6J. The peaks at the in- 
commensurate orientations do not appear in the angle dis- 
tributions at higher temperature, shown in Fig. 13(b). 

For graphite, the corrugation is only about 10 meV, and 
the low-temperature orbits significantly affect the friction 
only for temperatures below about 10 K. However, this 
effect can become important in a wide range of tempera- 
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tures for systems with larger potential corrugation, either 
due to higher normal load or to other substrates with 
hexagonal symmetry which exhibit larger corrugation, 
e.g. BN [20], which can be grown, for instance, on Rh [21]. 
In these cases, the frozen-orbit regime could provide an 
interesting new mechanism for low-friction sliding. 

In summary, we have investigated a simple model which 
nevertheless contains crucial elements of solid lubrication, 
in particular the rotational dynamics of the flakes. We 
have shown that quasi-superlubric behavior is achieved for 
the sliding of plates with disordered multi-domain struc- 
ture lubricated with graphite flakes. We have shown that 
rotational dynamics of the lubricant flakes, which destroy 
super lubric sliding of single flakes on a surface, are not 
necessarily detrimental in the actual geometry of solid lu- 
brication. 
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